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We study how the #-term is affected by interactions in certain one-dimensional gapped sys- 
tems that preserve charge-conjugation, parity, and time-reversal invariance. We exploit the 
relation between the chiral anomaly of a fermionic system and the classical shift symmetry 
of its bosonized dual. The vacuum expectation value of the dual boson is identified with 
the value of the #-term for the corresponding fermionic system. Two (related) examples 
illustrate the identification. We first consider the massive Luttinger liquid and find the 
#-term to be insensitive to the strength of the interaction. Next, we study the continuum 
limit of the Heisenberg XXZ spin-1/2 chain, perturbed by a second nearest- neighbor spin 
interaction. For a certain range of the XXZ anisotropy, we find that we can tune between 
two distinct sets of topological phases by varying the second nearest-neighbor coupling. 
In the first, we find the standard vacua at 9 = 0, n, while the second contains vacua that 
spontaneously break charge-conjugation and parity with fractional 9/% = 1/2,3/2. We 
also study quantized pumping in both examples following recent work. 
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1. Introduction 

1.1. Generalities 

Certain insulators can be distinguished by the presence or absence of a topological 
term in their low-energy effective actions [1]] . In even spacetime dimensions the topological 
term is a (9-term, while in odd dimensions, it is the Chern-Simons form. (By topological, 
we mean a term that is independent of the background metric.) More generally, non- 
interacting insulators are classified via the homotopy class of their gapped Hamiltonian 

HI- 

In this note, we are concerned with gapped systems defined in two spacetime dimen- 
sions and we study how the (9-term depends upon certain types of interactions. In the 
continuum, the canonical example of a theory whose low energy action contains a #-term 
is provided by a single massive Dirac fermion ift minimally coupled to a background U(l) 
gauge field A„. The Lagrangian has the form 



J d 2n x [i>(id^ - eA^ - m)^] , (1.1) 



where e, m are the charge and mass of the fermion, and 2n is the spacetime dimension. 
The matrices 7 M satisfy the algebra, {7^,7^} = 2?^, with rj^ = diag(l, —1, —1), and 
tp := tp^"f°. Call the field strength of the gauge field, F. Then, the 0-term, obtained in the 
low energy Lagrangian for the above example, is given by the quantity 

07(F) = - J d 2n x(eF/2n) n 1 (1.2) 

where F n := e il3x --- %n3n Fi x j 1 ■ ■ ■ Fi n j n . 9 takes a value equal to either zero or 7r depending 
upon the sign of the mass m. When 9 = the insulator is said to be trivial, while the 
insulator is called non-trivial when 9 = n. 

In general, 1(F) is integral when the theory is studied on a compact, oriented even- 
dimensional manifold.0 When the manifold is non-compact, quantization is a bit more 
subtle. In two dimensions, this condition is equivalent to the assumption that the back- 
ground U(l) gauge field is compact and the energy of a gauge field configuration is finite. 
In four dimensions, quantization of 1(F) follows from the assumption that the electro- 
magnetic field configuration is produced by electrically and magnetically charged point 



1 The factor of 1/2 ensures that the minimum value of 1(F) = 1 since we only consider systems 
on manifolds that admit a spin structure, i.e., fermions; otherwise, the 1/2 should be removed. 
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particles - Dirac quantization.ll Quantization of 1(F) implies that 9 has period 2tt since 
it is the quantity exp(idI(F)) that appears in the path integral. 

For a given spacetime dimension, 1(F) is odd under two of the discrete spacetime 
symmetries: charge-conjugation, parity, or time-reversal (C, V, T). Given that 9 has period 
27r and 1(F) is integral, the discrete symmetry under which 1(F) is odd folds the 2n interval 
in half; the two fixed points of the discrete symmetry are at 9 = 0, n. Thus, the positive 
and negative mass Dirac theories represent the two distinct classes of insulators invariant 
under the symmetries that change the sign of 1(F). 

To determine the discrete symmetries under which 1(F) is odd, it is necessary to 
know the transformation of the fields and coordinates. Charge conjugation does not act 
upon the coordinates, however, it changes the sign of all components of the gauge field, 
— > —A^. Parity is a reflection about the origin of all spatial directions Xi — >■ — Xi and 
similarly reverses the sign of all spatial components of the gauge field A{ — » —Ai. Time- 
reversal is an anti-unitary transformation, i — > —i, that transforms t — > — t and Ai — > —Ai. 
For example, in two and four spacetime dimensions, 1(F) transforms as (— , — , +) and 
(+,—,—), respectively, under (C,V,T), where ± refers to a quantity that is even/odd 
under the particular discrete symmetry. 

In perturbation theory, a non-zero #-term is found in the effective action for a back- 
ground gauge field coupled to a slowly varying charge density field after integrating out a 
massive Dirac fermion ||. The necessity of the charge density field should be contrasted 
with the situation in an odd number of spacetime dimensions. In this case, no charge 
density field is required and the Chern-Simons form is found at one-loop by integrating 
out a massive Dirac fermion. In odd dimensions, a component of the gauge field plays a 
role similar to that of the charge density field; this observation was used in the dimensional 
reduction technique of [|l| in deriving the existence of a #-term. A second approach for 
determining the #-term for massive fermionic systems is reviewed in the paragraph below. 
In Sec. 2, we show more explicitly how a #-term is generated in a related example. 

The following argument relies upon the U(l) chiral Adler-Bell-Jackiw anomaly ||. 
This argument is well known [[T],[7|,|8| and we review it for completeness. Denote the product 
of all gamma matrices, 72^+1 : = £ n 7o-"72«-i, and complexify the fermion mass term as, 

1 _ 1 
mipip = -mip(l + 72n+l)'^/ , + t^ 7 ™*^ 1 ~ ^2n+i)^- (1.3) 

2 For a recent discussion of the quantization of 1(F) in 3+1 dimension, please see [||. 
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Under a chiral transformation, tp — > exp(zcry2n+i)' i /'; the complex mass parameter rotates 
m — > exp(2za)m. Because the path integration measure is not invariant under the chiral 
rotation, the quantum action also shifts by a term proportional to 2aI(F) ||. (That the 
action for the Dirac fermion shifts by a #-term with 9 = 2a is exact at one-loop.) However, 
physical quantities cannot depend upon the choice of variables in the path integral. Ob- 
servables can only depend upon the product M = mexp(— z#) with 9 being the coefficient 
of an 1(F) that may be present in the tree-level action. Thus, the theory ( |1 . 1|) with posi- 
tive mass is different from the one with negative mass with vanishing tree-level 9 in both 
theories because M differs by a sign. However, it is possible to rewrite the Lagrangian 
for the negative mass theory by noting — m = mexp(z7r), and so the two massive Dirac 
Lagrangians only differ by a term equal to ttI(F). Once the massive fermion is integrated 
out, the effective Lagrangians for the background gauge field will differ by the presence of 
the term ttI(F). 

Please see || for a very nice discussion about how anomalies in various dimensions 
can be used to classify gapped systems. 

1.2. Motivation and Outline 

It is interesting to generalize these ideas to interacting topological phases [p70| , p7T| , |T2|| . 
In the context of the massive Dirac theory, a generic interaction does not transform simply 
under a chiral transformation and so the above argument does not apply in general. In 
particular, a chiral rotation can change the sign of an interaction term and the low energy 
properties of a theory can be sensitive to the sign change. However, even if the interaction 
is invariant under a chiral transformation, the change of measure of the path integral, as 
inferred from an anomalous axial current, can naively depend upon the interaction. As 
an example of the latter phenomenon, in Sec. 2 we study the 1+1 dimensional massive 
Luttinger liquid; here, the axial anomaly equation explicitly depends upon the interaction 
strength and the Fujikawa analysis is complicated by the presence of the marginal current- 
current interaction. In general, however, a relevant interaction introduces a new mass scale 
and may dramatically affect the low energy physics. 

There are two approaches towards the study of interacting topological insulators. 
The first is to consider a massive theory and add arbitrary interactions that maintain the 
symmetries of interest. We do this in Sec. 2. The second approach is less direct and 
generally rather difficult to study in practice. Begin with a classically gapless Lagrangian 
that admits a massive perturbation and ask whether or not the particular deformation 
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results in a physically interesting topological phase. We take this approach is Sec. 3. 
There need not be complete overlap with the various phases obtained (theoretically) via 
the above two approaches. 

Because the models we discuss and the techniques by which we study them are well 
known, let us summarize what is new. Using bosonization to study interacting 1+1 di- 
mensional fermionic systems, we identify the vacuum expectation value of the dual boson 
with the value of the (9-term. Bosonization allows us to determine the #-term directly. 
We familiarize ourselves with this identification in the context of the massive Luttinger 
model. Next, we relate the dimer and Neel phases of the Heisenberg XXZ spin-1/2 chain 
perturbed by a second nearest-neighbor interaction to the value of a corresponding #-term: 
the dimer and Neel phases contain the 9 = n/2,3n/2 and 9 = 0, n vacua, respectively. 



Finally, we study adiabatic charge pumping |T3,p^[T5| in both models. The novelty in the 



XXZ chain lies in the fact that it is necessary for the varied parameter to wind twice about 
the origin in order to transfer a full unit of charge across the system; a single winding is 
analogous to varying between a positive and negative mass theory in the context of the 
Luttinger liquid. 

This paper is organized as follows. In Sec. 2, we study the massive Luttinger liquid 
as an example of a 1+1 dimensional interacting topological phase. We find that the Z2 
classification by the #-term is independent of the strength of the interaction. In Sec. 3, 
we discuss the continuum limit of the Heisenberg XXZ spin-1/2 chain, perturbed by a 
second nearest-neighbor interaction. For a certain range of the XXZ anisotropy, we find 
that we can tune between two distinct sets of topological phases by varying the second 
nearest-neighbor coupling. In the first set, we find the standard 9 = 0, n vacua, while the 
second set contains vacua that spontaneously break charge-conjugation and parity with 
fractional 9 = n/2,3n/2. Adiabatic variation of certain parameters in both models results 
in charge transport. In contrast to the massive Luttinger liquid |l3|Jl5[| , the XXZ chain 



transitions to the 9 = n vacuum after a single closed orbit in parameter space and only 
results in an half- unit of charge transfer. Two windings are necessary in order to return 
to the original vacuum and so the total charge transferred is integral. We summarize in 
Sec. 4. 



2. Massive Luttinger Liquid 

2.1. Invariance of 9 

A simple example where the interactions can be treated exactly is provided by the 



massive Luttinger liquid fll(||17| , |18| . The Lagrangian 

L = J d\^^7 M -e^ 7 ^-^(ma^ (2.1) 

describes a 1+1 dimensional spinless fermion ijj of charge e coupled to a background 
gauge field A^. The constant charge density field m couples the left and right handed 
components of ip and functions as the mass parameter. The Lagrangian breaks both C and 
V when m has a non-zero imaginary part. The current-current interaction is parameterized 
by the coupling g. A more general quartic interaction will be considered in Sec. 3. 

At finite g, it is simplest to study (|2.1| ) after bosonization of the fermion [fl8 |; we 
introduce a scalar field 6 that satisfies 



l + <7/7r = 47r//3 2 , 
In 



(2.2) 



+ 73)^ = ^R^L 



2(3< 



exp(^0), 



where the magnitude of the coefficient of the exponential in the final equality in (|2.2|) is 
chosen for convenience. Substituting ( |2.2|) into (|2.1[) , and after one integration by parts in 
the last term, we obtain 



d 2 x 



1 



(9, 



711 CRA\ m 

exp{ip<p) — 



2(3 



2/3' 



exp(-i/30) + ^.^d„A v 
Ztv 



(2.3) 



The integration by parts in the last term of ( |2.3| ) is required if the action is to be gauge- 
invariant on a space with boundary. Further, this form of the action correctly displays 
the chiral anomaly associated with shifts of the <p field or, in the fermionic language, the 
Jacobian of the path integration measure under a chiral transformation. More precisely, a 
chiral rotation of if) — > exp(ia~f3)i() is equivalent to a shift of <p by a and the last term in 
( |2 . 1| ) correctly reproduces the shift of the action by 1(F) under this transformation. We 
also note that </> — > —<fi + n under a C transformation, while <p changes sign under a V 
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transformation, and is neutral under T . (Although C effectively takes m — > —m, it also 
shifts 9, thus preserving mexp(-iO).) 

Expanding the exponential exp(±z/3</>), (3 2 acts as a coupling constant for the leading 
quartic scalar interaction, assuming for the moment that m is real. Thus, the duality 
( p.2|) interchanges a strongly coupled fermion theory for a weakly coupled boson theory 
and vice versa. When (3 2 = 4tc, g vanishes and ( |2.3|) relates a free massive fermion to an 
interacting scalar field. Through its definition in terms of /3, the interaction strength g of 
the four-fermion interaction controls the radius of the <p boson; <p is periodically identified 
under — >• + 2 n/ft. Invariance of exp(iL) under a shift of <p by its period is consistent 
with the integrality of 1(F); in general, however, integrality of 1(F) is distinct from the 
requirement that exp(iL) be invariant under periodic shifts of (p. For the remainder, we 
restrict (p to the interval [0, 2tt/(3). 

The scaling dimension of exp(i(3(p) equals (3 2 /4tc or (1 + g/n) -1 in terms of the 
fermionic interaction strength g. Attractive interactions between the fermions increase 
the dimension of the exponential operator, while repulsive interactions have the opposite 
effect. (3 2 > 8n or g < —n/2 parameterize a line of fixed points that terminates when 
/3 2 = 8tt at the SU(2) 1 WZW critical point |9f. Thus, we restrict to the regime where 
exp(i(3<p) is relevant. 

It is convenient to absorb the phase of the mass m by redefining (p. Writing m = 
\m\ exp(ia), we rename <p = (p + a/ (3. With this redefinition, the Lagrangian becomes 

L = J d 2 x\^{d^) 2 - M C os(/30) + ^(0-a//3)e^a M ^]. (2.4) 

At energies low compared to \m\ the scalar field is frozen and minimizes its potential energy 
at (<p) = n/(3 or, in the original coordinate, (<p) = (n — a)/ (3. In the low energy effective 
Lagrangian 

Lev = ^-y^ J d 2 x(e e^F^/2n), (2.5) 

the dependence upon (3 has completely disappeared; the effective Lagrangian is indepen- 
dent of the strength of the current-current interaction. We see explicitly that the low 
energy Lagrangian equals (n — a)I(F). The C and V symmetric points lie at a = 0,7r; 
namely, when m > or m < O.li When a differs from these two values, the theory breaks 
C and V. These results are independent of the interaction strength. 



3 We remark that identifying the ±m theories with the 9 = 0, ir vacua is a convention. An 
invariant statement requires a specification of a tree-level #-term - tuned to zero in this example - 
as discussed in the Introduction. In the absence of this, we can only say that the effective actions 
for the two theories differ by a 6 = it term. 
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2.2. Adiabatic Variations 



In this section, we re-derive the results recently obtained in |jT5|l , and discussed previ- 
ously in slightly more generality in [|14| ,|T| , concerning the Thouless charge pump |l3j . This 
review provides intuition for a similar study in the context of the Heisenberg XXZ spin-1/2 
chain discussed in Sec. 3. As an aside, we mention how the chiral anomaly equation is 
related to the conductance of a one-dimensional wire. 



It is natural to consider two currents constructed from the fields in (|2.1|) : (1) the vector 
current Jy = ^^if) and (2) the axial current = '^/>737' x '^/ , '. Their integrated charge 
densities describe, respectively, the sum and difference of the number of left and right 
moving fermions. If both currents were conserved, the number of left and right moving 
fermions would be separately conserved. However, the axial current is not conserved 
because of the fermion mass term and the minimal coupling to the background gauge field. 
The bosonization rules ( |2.2| ) allow us to rewrite these currents in terms of the scalar 
Jy = j-e^ v d u (f) while = J^d^rf*. The vector current is identically conserved while the 
axial current, given its definition in terms of <p, obeys an equation that follows directly 
from the (p equation of motion, 



d,J% = ^1 sin(^) + 'erd^. (2.6) 



\ m \ ■ i n~,\ /3 2 e 

In the fermion language, the first term on the righthand side of ( |2.6| ) is the explicit violation 
of the axial symmetry by the mass term while the second term obtains from a variation of 
the fermion path integration measure. Note that the coefficient of e^ v d^A v depends upon 
the interaction through /3 2 /47r. 

Consider first an adiabatic change of sign of the fermion mass when the background 
field is set to zero. In particular, consider the arc parameterized at non-zero \m\ by 
adiabatically varying the phase a with respect to time. In the vacuum, the vector current 
{Jy) = 2^^ V {dv4>) = ~^ e ^ u dv a - Therefore, the amount of charge that passes through 
an arbitrary point on the line, 

e f°° , „ e 



oo 



AQ = e altJy = — alt d t a = —Aa. (2.7) 

2tt J_ nn 2tv 



Choosing Aa = 27r, i.e., a loop about the origin, ( |2.7[ ) is simply a restatement of Thou- 



less' result on quantized charge transport in one spatial dimension [13] generalized to the 



Luttinger liquid as nicely demonstrated recently in |jl5| |. Adiabatic variation of the phase 
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of the mass term about the origin results in transport of a single unit of charge (modulo 
one) through the system. If the loop does not enclose the origin, Aa = and no charge is 
transported. ( |2.7| ) gives an explicit example with interactions of the phenomenon discussed 
in [|| where two distinct topological insulators are connected via a path in parameter space 
that preserves a non-zero excitation gap, but does not necessarily maintain the discrete 
symmetries of the two classes of insulator. 

Contact with the work of pO| , pT| ,|5| can be made by considering instead an interpolation 
in space instead of time between two values of m. It is sufficient to consider an interpolation 
of the phase oof m. The charge carried by a soliton, 

Q = e [ dxJ^ = e [ dx—d x (4>) = — Aa. (2.8) 
J J 2n 2h 

If the soliton interpolates between positive and negative values of the mass parameter, 
then Aa = n and so Q = e/2 [p0 |. However, in general, the change of the phase of m, and 
therefore the corresponding charge carried by the soliton, can be arbitrary ||. 

This provides an heuristic way to understand why the charge transported for any 
closed loop in m space is integral. Consider a finite length chain and imagine introducing 
a soliton or kink near the left end that interpolates between the two vacua at ±m. This 
interpolation is achieved by varying the phase of m as above in space as opposed to time. 
The charge carried by the soliton, Q = e/2. Closing the loop by returning to the positive 
m— axis results in the transport of a second e/2 soliton. If we had instead reversed the path 
in m space so that the origin was not enclosed, it would have been equivalent to taking a 
kink from the left to the right end and then back again with zero net charge transfer. 

As an aside, we mention a result that follows from the axial current equation ( |2.6|) . 
Set m = and allow the external electric field & v d^A v to vary adiabatically with time. 
In particular, adiabatically vary A x while keeping At = 0. The system is placed on a 
spatial circle of length L and ip is taken to be periodic around L (more general boundary 
conditions are also possible). A general constant background A x cannot be removed by 
a gauge transformation as this would ruin the assumed periodicity of ijj since a gauge 
transformation multiplies the fermion field by exp(iece) with a the gauge transformation 
parameter. However, when A x = 2nn/eL, for n an integer, it is possible to perform a 
gauge transformation with the choice a = A x x to remove the background field. In other 
words, a gauge field configuration A x = 2T\njeL is physically equivalent to a configuration 
without any background field. 
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Let A x adiabatically vary from zero to 2tt JeL so that the system returns to its orig- 
inal ground state after the variation with, perhaps, a slight relabeling of its spectrum. 



Integrating (|2.6| ) with m = over spacetime, we find 

?2„ r i&. 



N L -N R = J dtdx d^ A = £| J dtdx d t A x = ^^AA X . (2.9) 



Because AA X = 27r/eL, the difference Nl — Nr changes by Note that if the interaction 



2tt' 

?2 



g is tuned to zero, /3 = 4% and the difference Nl — Nr = 2 reflecting the fact that one 
left-handed fermion has emerged from the fermi sea, while one right-handed fermion has 
been captured. Such a flow in the spectrum ensures the charge of the vacuum remains 
the same after the adiabatic variation. However, a non-zero current obtains after the 
variation since there is an imbalance between the number of left and right movers. For 
arbitrary interaction, that the difference is proportional to (3 2 /4n is a manifestation that 

o2 2 

the conductance of an interacting 1+1 dimensional chain equals 4^:7^ (assuming Luttinger 
liquid leads) 



3. Generic Interactions 

In this section, we study slightly more generic interactions that preserve C and V. We 
focus on a one- dimensional system that is gapped in the IR, not because of a tree-level mass 
term, but because of interactions.^ We view this as an alternative approach to building 
insulators with interesting topological properties. 

3.1. 6 -Vacua 

The motivation for the model considered in this section is the Heisenberg XXZ spin- 
1/2 chain described by the microscopic Hamiltonian, 

H = 3 E [ S * • S ^ + S " • #+1 + 9' St ■ S? +1 ] , (3.1) 



1 



with J, g' > 0. The spin operators obey the SU(2) algebra, [Sf,Sj] = ie abc Sj for spins 
on the same lattice site, j. g' parameterizes a line of critical points for < g' < 1. The 
system enters a massive phase when g' > 0. (The system is ferromagnetic when g' < 0.) 



See p5| | for a related discussion in 3+1 dimensions where the interactions are due to disorder. 
In general, disorder may be treated exactly in 1+1 dimensions [p6[|. 
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It is well known that the continuum limit is described by the Lagrangian (for a review, see 



d 2 x 



^{id^ ~ - |((^) 2 " 2((^l) 2 + (^R?) ■ (3.2) 



where the subscripts L, R refer to the left and right moving parts of the fermion. In obtain- 
ing ( p.2| ) , we have dropped unimportant constants and absorbed any renormalization of the 



Fermi velocity by a redefinition. Note that we have introduced in (|3.2|) a background gauge 
field minimally coupled to the number current of the Jordan- Wigner fermions. Fermion 
number is identified with the S z eigenvalue in the spin language. The interaction terms 
are parameterized by the coupling g which is linearly related to the XXZ anisotropy pa- 
rameter, g' . In ( |3.2|) , the first interaction term is simply the current- current interaction 
whose effects we previously studied in Sec. 2. The second interaction term, naively zero by 
the Pauli principle, is responsible for dynamical mass generation. The product should be 
understood as the leading term in the operator product expansion (OPE) of the operator 
(ipx L ipL,R) with itself. This second term is the spin-less analog of Umklapp scattering 
|2^,[29],|3(| and we refer to it as the Umklapp term. 

Bosonization proceeds as in Sec. 2 by introducing a scalar field <fi satisfying (|2~2j). The 
resulting Lagrangian 



/ [l^ 2 + J? cos (W) + If^M*] • ( 3 - 3 ) 



We assume that we are in the regime (3 2 < 2n where the Umklapp term is relevant. (The 
precise value in terms of microscopic parameter g' of the XXZ chain at which this transition 
occurs can be obtained and supports the qualitative prediction of the sine-Gordon model 
that a transition to a massive phase occurs for some finite g ~ 7T.) The low energy limit of 
( jOl) is achieved by taking g — >■ oo. The Umklapp potential freezes the scalar field at one 
of the two vacua, (</>) = 0, 7r//3. These two vacua are Z2 chiral partners and correspond to 
9 = 0, 7r, respectively. 

Notice that in contrast to the example in Sec. 2 where the tree level mass term only 
allowed a single vacuum for a particular sign of the cosine interaction (assuming the range 
of (p was restricted to the interval [0, 27r//3)), the two topologically distinct vacua for this 
model both occur for one sign of the cosine interaction. Of course, this is due to the doubled 
periodicity of the cosine potential. Note there do not exist additional non-trivial vacua 
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when g < 0. The reason is that the dimension of the Umklapp operator is determined by 
g and so when g < tt, cos(2(3(p) is irrelevant and decouples from the low energy physics. 

It is possible to find additional vacua when the coefficient of the cosine is negative, 
however, it is necessary to independently vary the dimension of the Umklapp operator and 
its coefficient. In terms of the motivating XXZ chain, this is achieved by adding a second 



nearest-neighbor spin interaction [31 



5H = A^S,-S l+2 . (3.4) 

i 

The bosonic version of the continuum Lagrangian 

L = j [l(d^) 2 + cos(2/3</>) + ^-^d»A v \ , (3.5) 

where (5 2 /Att = (1 + (g + 2A)/7r) _1 . Choose g > it so that the Umklapp operator is 
relevant even when the second nearest-neighbor coupling vanishes. The massless Luttinger 
liquid separates at A = g/6 two distinct sets of massive phases . When A < g/6, the 
vacua lie at (</>) = 0,7r//3, while when A > g/6, (</>) = n/2/3, 3tt/2(3 and the corresponding 
9 = n/2,3n/2. The two vacua in the regime A > g/6 spontaneously break C and V. This 
follows from the fact that 1(F) is integral. 

If 1(F) were even, for example, the range of 9 would be reduced to the interval [0, n) 
and the two vacua at 9 = 7r, 37r/2 would be identified with their tt / (3 partners and effectively 
removed. The remaining vacua would become the two C, V invariant vacua - the vacuum 
with non-zero vacuum expectation value having fractional 9 = n/2. We are not aware 
of an argument that would require 1(F) to be even. Note that invariance of the total 
action with respect to a shift of <p that preserves the cosine interaction is equivalent to the 
requirement that 1(F) be even, however, quantization of 1(F) and the shift symmetry of 
a potential need not be related. 

In terms of the spin-chain, the transition at A = g/6 separates two distinct sets of 
phases. The vacua that preserve C, V correspond to the Neel phase while the dimerized 
phase occurs for A > g/6, given a sufficiently large g > it. 

Finally, we remark that ( |3.5| ) contains the most relevant interactions consistent with 
the symmetries of the spin chain. One possible operator of lower dimension is cos (/?</>), 
obtained from an alternating interaction in the spin model, but it can be ignored since it 
violates the discrete Z2 lattice symmetry, <fi — > <fi/(3. There is one other possibly dangerous 
operator, namely, cos(/3</>) where (f> is defined as follows. Expand <fi = <pL + 4>R l e ft 
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and right moving modes. We define (p := 4>l — 4>r- (It can be identified with the T-dual 



of <f>.) It is possible to ignore cos(/3</>) as long as the spin-chain Hamiltonian (3.1) enjoys 
the U(l) rotation symmetry about the z-axis since this acts as <fi — > <fi + a where a is an 
arbitrary constant. 

3.2. Adiabatic Variation 

Following the method discussed in Sec. 2, it is now possible to examine the nature of 
charge transport under adiabatic variation of parameters in the XXZ spin-chain. Just as it 
was necessary to allow the mass parameter in Sec. 2 to be complex, we take the coefficient 
of the cosine interaction discussed in this section to be complex. This is done as follows. 

The continuum limit of the XXZ spin-chain is parameterized by two couplings, g, A 
which we allow to be complex, g + 2 A determines the dimension of the cosine interaction 
while and g — 6A is the coupling constant of the interaction. Thus, we must keep the linear 
combination g + 2A real while allowing g — 6A to take complex values. Hermiticity of the 
Lagrangian is maintained by adding in the Hermitian conjugate interaction. Writing g = 
91+W2 and A = A1-MA2, for Oj, \ real, we consider arbitrary values of g\ — 6A1 +i(g 2 — 6A2) 
while requiring g 2 = — 2A2 so that the dimension of the operator is real. The relation 
between g 2l X 2 reduces the dimension of the parameter space by one. We take gi,g%, Ai to 
parameterize the continuum limit of the XXZ chain. 

In order to simplify the presentation, let us define, g = g\ + 2Ai and pexp(ia) = 
g 1 - 6A1 + 4ig 2 , with < p = y/( gi - 6A1) 2 + 16#| and a = tan- 1 (4^ 2 /(^i - 6A1)). The 
resulting bosonic Lagrangian 

L= [ d 2 x\l(d^) 2 + pcos(2p ( j ) + a) + ^ ( j ) e^d^AX (3.6) 
J 12 Zn J 

where Air / (3 2 = 1 + g/%. At low energies or p — > 00, the vacua lie at (</>) = —a/2/? + nir/j3, 
for integer n. 

We consider the path obtained at non-zero p by varying a : — > 2n. Following Sec. 2, 
the amount of charge that is transferred across the system 

AQ = -e dt£-d t ((t>) = — dtd t a= —Aa. (3.7) 
J -00 27r 2tt J_ 00 4tt 



The factor of 1/2 on the right hand side of ( |3.7| ) is due to the doubled periodicity of the 



cosine interaction. When Aa = 2tt, an half-integral charge is transferred across the system. 
The system, however, has not returned to its initial state, but has transitioned to its Z 2 
chiral partner. It is necessary to wind twice about the origin, a : — > 4/T, in order to 
return to the original vacuum. Doing so, we find integral charge transport. 
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4. Discusion 

We considered the effects of certain interactions on the #-term for one-dimensional 
gapped systems. Our examples were simple and well known. What is new is our viewpoint 
on the transitions described by tuning of various parameters in the models. We made use 
of abelian bosonization techniques to determine the value of the #-term in the low energy 
action. In particular, the semi-classical minimum of the potential for the dual boson was 
shown to directly determine the value of the (9-term. 

The first model we studied was that of the massive Luttinger liquid. When the 
interactions are taken to vanish, the model describes free massive fermions. The sign of 
the mass determines the topological phase in which the system lies. The two phases are 
distinguished by the value of a 6>-term. We found these two values of 6 to be stable to 
arbitrary current-current interactions. 

In the course of studying this model, we made contact with the recent work |15|j . By 
adiabatically varying the phase of the mass term, a unit of charge is pumped through the 
system. In the massless limit, we also noticed that the axial anomaly equation provided a 
simple way to infer the conductance of a one-dimensional interacting wire. 

The second model we studied was inspired by the Heisenberg XXZ spin-1/2 chain 
perturbed by a second nearest-neighbor interaction. For vanishing second nearest-neighbor 
interaction, we found two vacua related to each other by a Z<i chiral translation. Phases 
with fractional 6 /tc = 1/2, 3/2 obtained when this interaction was turned on with sufficient 
strength, however, these vacua spontaneously broke both C and V. In XXZ chain, we 
studied adiabatic charge pumping and found, in contrast to the model studied in Sec. 2, 
it was necessary to perform two cycles in order to return to the same vacuum at the end 
of the variation. Again, an integral unit of charge was pumped during the variation. 

Part of the initial motivation in studying simple one-dimensional systems amenable to 
bosonization techniques was to have a better understanding of the (9-term in the context 
of various topological phases of matter. In particular, we were interested in studying 
interactions in the controlled environs of one-dimensional systems. This work provides 
context for further exploration of phases of matter where we can hope to find fractional 
values of 9/n preserving C, V, T [P2]|7|. 



13 



Acknowledgements 

It is a pleasure to thank F. Benini, F. Burnell, K.-T. Chen, C. Herzog, C. Hoyos- 
Badajoz, S. Kachru, P. Lee, J. McGreevy, A. Scardicchio, and B. Swingle for discussions 
and encouragement. We acknowledge the hospitality of the Stanford Institute for The- 
oretical Physics, the Galileo Galilei Institute for Theoretical Physics, and Oxford Uni- 
versity while this work was in progress. We were supported in part by funds provided 
by the U.S. Department of Energy (D.O.E.) under cooperative research agreement DE- 
FG0205ER41360. 



14 



References 



[1] Xiao-Liang Qi, Taylor L. Hughes, and Shou-Cheng Zhang, "Topological field theory 
of time-reversal invariant insulators," Phys. Rev. B 78, 195424 (2008) 

[2] Alexei Kitaev, "Periodic table for topological insulators and superconductors," 
^rXiv:0901 .26561 [cond-mat]. 

[3] Shinsei Ryu, Andreas Schnyder, Akira Furusaki, Andreas Ludwig, "Topological insu- 
lators and superconductors: ten-fold way and dimensional hierarchy," New J. Phys. 
12, 065010 (2010). 

[4] M.M. Vazifeh, M. Franz, "Quantization and 27r Periodicity of the Axion Action in 

Topological Insulators," |arXiv:1006.3355| [cond-mat]. 
[5] Jeffrey Goldstone and Frank Wilczek, "Fractional Quantum Numbers on Solitons," 

Phys. Rev. Lett. 47, 986 (1981). 
[6] Stephen L. Adler, "" Axial- Vector Vertex in Spinor Electrodynamics," Phys. Rev. 177, 

2426 (1969); J. S. Bell and R. Jackiw, ""A PCAC puzzle: r)° 77 in the a-model" 

Nuovo Cim.A60:47-61,1969. 
[7] Carlos Hoyos-Badajoz, Kristan Jensen, Andreas Karch, "A Holographic Fractional 

Topological Insulator," |arXiv:1007.3253| [hep-th]; Andreas Karch, Joseph Maciejko, 

Tadashi Takayanagi, "Holographic fractional topological insulators in 2+1 and 1+1 

dimensions," |arXiv: 1009 .29911 [hep-th] 
[8] Shinsei Ryu, Joel E. Moore, Andreas W. W. Ludwig, "Electromagnetic and grav- 
itational responses and anomalies in topological insulators and superconductors," 

|arXiv: 1010.09361 [cond-mat]. 
[9] Kazuo Fujikawa, "Path Integral for Gauge Theories with Fermions," Phys. Rev. D 21, 

2848,(1980). 

[10] Lukasz Fidkowski and Alexei Kitaev, "Effects of interactions on the topological clas- 
sification of free fermion systems," Phys. Rev. B 81, 134509 (2010); Lukasz Fidkowski 
and Alexei Kitaev, "Topological phases of fermions in one dimension," |arXiv: 1008.4138 
[cond-mat] . 

[11] Xie Chen, Zheng-Cheng Gu, Xiao-Gang Wen, "Classification of Gapped Symmetric 
Phases in ID Spin Systems," |arXiv:1008.3745| [cond-mat]; Zheng-Cheng Gu, Zhenghan 
Wang, Xiao-Gang Wen, "A classification of 2D fermionic and bosonic topological 
orders," |arXiv:1010.1517j [cond-mat]. 

[12] Ari M. Turner, Frank Pollmann, Erez Berg, "Topological Phases of One-Dimensional 
Fermions: An Entanglement Point of View," |arXiv: 1008.4346| [cond-mat]. 

[13] D. J. Thouless, "Quantization of particle transport," Phys. Rev. B 27, 6083 (1983). 

[14] Liang Fu, C.L. Kane, "Time Reversal Polarization and a Z2 Adiabatic Spin Pump," 
Phys. Rev. B 74, 195312 (2006). 



15 



[15] Erez Berg, Michael Levin, and Ehud Altman, "Quantized pumping and phase diagram 

topology of interacting bosons," |arXiv:1008.1590| [cond-mat]. 
[16] Sin-itiro Tomonaga, "Remarks on Bloch's Method of Sound Waves applied to Many- 

Fermion Problems," Prog. Theor. Phys. Vol. 5 No. 4 (1950). 
[17] J. M. Luttinger, "An Exactly Soluble Model of a Many?Fermion System," J. Math. 

Phys. 4, 1154 (1963). 

[18] S. R. Coleman, "Quantum sine-Gordon equation as the massive Thirring model," 

Phys. Rev. D 11, 2088 (1975). 
[19] Daniel J. Amit, Yadin Y. Goldschmidt, and G. Grinstein, "Renormalisation group 

analysis of the phase transition in the 2D Coulomb gas, Sine-Gordon theory and XY- 

model," J. Phys. A: Math. Gen. 13 585 (1980). 
[20] R. Jackiw and C. Rebbi, "Solitons with fermion number 1/2," Phys. Rev. D 13, 3398 

(1976). 

[21] A. J. Heeger, S. Kivelson, J. R. Schrieffer, W. -P. Su, "Solitons in conducting poly- 
mers," Rev. Mod. Phys. 60, 781850 (1988). 

[22] W. Apel and T. M. Rice, "Combined effect of disorder and interaction on the conduc- 
tance of a one-dimensional fermion system," Phys. Rev. B 26, 7063-7065 (1982). 

[23] Dmitri L. Maslov and Michael Stone, "Landauer conductance of Luttinger liquids with 
leads," Phys. Rev. B 52 5539 (1995). 

[24] Ronny Thomale and Alexander Seidel, "A minimal model of quantized conductance 
in interacting ballistic quantum wires," |arXiv:1005.4228| [cond-mat]. 

[25] H.-M. Guo, G. Rosenberg, G. Refael, M. Franz, "Topological Anderson Insulator in 
Three Dimensions," |arXiv:1006.2777| [cond-mat]. 

[26] T. Giamarchi and H. J. Schulz, "Anderson localization and interactions in one- 
dimensional metals," Phys. Rev. B 37, 325 (1988). 

[27] E. Fradkin, Field Theories of Condensed Matter Systems" (Addison- Wesley, Redwood 
City, 1991); I. Affleck, "Field Theory Methods and Quantum Critical Phenomena", 
Fields, Strings and Critical Phenomena, p. 563-640, (ed. E. Brezin and J. Zinn- Justin 
North- Holland, Amsterdam, 1990), proceedings of Les Houches Summer School, 1988, 
invited lectures. 

[28] F. D. M. Haldane, "General Relation of Correlation Exponents and Spectral Properties 

of One-Dimensional Fermi Systems: Application to the Anisotropic S=l/2 Heisenberg 

Chain," Phys. Rev. Lett. 45, 1358 (1980). 
[29] J. L. Black and V. J. Emery, "Critical properties of two-dimensional models," Phys. 

Rev. B 23, 429 (1981). 
[30] M. P. M. den Nijs, "Derivation of extended scaling relations between critical exponents 

in two-dimensional models from the one-dimensional Luttinger model," Phys. Rev. B 

23, 6111 (1981). 



16 



[31] F. D. M. Haldane, "Spontaneous dimerization in the S=l/2 Heisenberg antiferromag- 
netic chain with competing interactions," Phys. Rev. B 25, 4925 (1982). 

[32] Joseph Maciejko, Xiao-Liang Qi, Andreas Karch, Shou-Cheng Zhang, "Fractional 
topological insulators in three dimensions," |arXiv:1004.3628| [cond-mat]; Brian Swingle, 
Maissam Barkeshli, John McGreevy, T. Senthil "Correlated Topological Insulators and 
the Fractional Magnetoelectric Effect," arXiv:lUU5.1U76| [cond-mat]. 



17 



